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ABSTRACT 

In view of the recent measurement of the scintillation velocity for PSR J0737-3039, we examine the 
complete set of constraints imposed on the pulsar B natal kicks (magnitude and orientation) and 
predict the most favorable pulsar kick velocity and spin tilt for both isotropic and polar kicks. Our 
analysis takes into account both currently unknown parameters: the orientation of the orbital plane 
on the sky (f2) and the radial component of the systemic velocity (V r ). Assuming that the system's 
peculiar velocity is entirely due to the second supernova explosion, we find that the system may have 
crossed the Galactic plane multiple times since the birth of the second neutron star and that the 
post-supernova peculiar velocity could have been as high as ~ 1200kms _1 . We also confirm the 
absolute lower and upper limits on the physical parameters derived in our earlier study. For specific 
combinations of the two unknown parameters f2 and V r , however, we find much tighter constraints on 
the pre-supernova binary configuration and natal kicks imparted to pulsar B, as well as on the age of 
system. Once fi is measured in the coming year, it will be straightforward to use the results presented 
here to further constrain the natal kicks and the spin-tilt predictions. We complete our comprehensive 
study and derive similar constraints and spin-tilt predictions for PSRB1534+12, where the only free 
parameter is V r . Lastly, for PSRB1913+16, we update the progenitor and kick constraints using 
the measured pulsar spin tilt and allowing for Roche-lobe overflow from the progenitor of the pulsar 
companion. 

Subject headings: Stars: Binaries: Close, Stars: Pulsars: General, Stars: Neutron, Stars: Pulsars: 
Individual: (PSR J0737-3039, PSRB1534+12, PSRB1913+16) 



1. INTRODUCTION 

The recent discovery of a unique, strongly relativis- 
tic binary pulsar (Burgay et al. 2003) that is also the 
first eclipsing, double pulsar system found in our Galaxy 
(Lyne et al. 2004) promises to provide us with more 
impressive tests of general relativity than ever before 
(Kramer 2004), to greatly improve our understanding of 
magnetospheric and pulsar emission physics (Lyutikov 
2004; McLaughlin et al. 2004; Demorest et al. 2004; 
Granot & Meszaros 2004; Arons et al. 2004), and possi- 
bly indicates a realistic chance for the detection of grav- 
itational waves from double-neutron-star (DNS) inspiral 
with the first-generation of ground-based interferometers 
(Kalogera et al. 2004; Kim et al. 2004). 

Apart from the binary pulsar and eclipsing behavior, 
PSR 0737-3039 sets the current records for DNS prop- 
erties: for example, it has the shortest spin and orbital 
period, smallest eccentricity, lowest magnetic field, and 
fastest apsidal motion. The measured current orbital 
parameters and masses along with an estimate of the 
system's age have allowed the careful account of its re- 
cent evolutionary history, since the time just before the 
supernova (SN) that formed pulsar B, the current com- 
panion to the millisecond pulsar A (Willems & Kalogera 
2004; Dewi & van den Heuvel 2004). In Willems & 
Kalogera (2004) (hereafter Paper I), we found that the 
pre-SN binary 1 most likely consisted of the first-formed 

Electronic address: b-willems, vicky@northwcstern.edu, 
m-henninger@alumni . northwestern.edu 

1 Throughout the paper, we use the terms pre- and post-SN to 
refer to the instants just before and right after the SN explosion 
forming the second NS. 



NS and a Roche-lobe-filling helium star with mass be- 
tween 2.1 Mq and 4.7 M & (cf. Dewi & van den Heuvel 
2004). The pre-SN orbital separation can be constrained 
between 1.36 Rq and 1.72 Rq, and the magnitude of the 
kick velocity between 60kms~ 1 and 1560 km s -1 . 

Most recently, scintillation observations of pulsar A 
have allowed the measurement of the two systemic ve- 
locity components in the reference frame of the orbital 
plane projected on the sky (Ransom et al. 2004). This 
measurement gives us the first glimpse of the kinematic 
evolutionary history (in addition to that related to bi- 
nary evolution). Therefore we extend the analysis pre- 
sented in Paper I and take a comprehensive approach 
to constraining the formation of PSR J0737-3039. The 
concept is similar to that adopted by Wex, Kalogera, & 
Kramer (2000) with two main differences: (i) at present 
there is no constraint on pulsar A's spin tilt from timing 
or polarization observations; (ii) the current kinematic 
measurements involve two free parameters, the orienta- 
tion of the orbital plane on the sky (fi) and the radial 
component of the systemic velocity (V r ). In 
this paper, we account for the complete set of equations 
describing the formation of PSR J0737-3039 and derive 
constraints on the progenitor of pulsar B and on the 
kicks (both isotropic and polar) imparted to it at birth. 
We also derive predictions for the possible orientation 
of pulsar A's spin axis with respect to the current or- 
bital angular momentum axis. We carefully analyze and 
present the dependence of these results on and V r and 
on the system's age. The measurement of f2 along with 
proper motion components is anticipated in the coming 
year. We show that such measurements can be used with 
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TABLE 1 

Parameters of PSR J0737-3039, PSRB1534+12, and PSRB1913+16. 



Parameter 


Notation 


PSR J0737-3039 a 


PSRB1534+12 b 


PSRB1913±16 C 


Right Ascension (J2000) 


a 


07 h 37 m 51.25 s 


15 h 37 m 09.96 s 


19 h 15 m 28.00 s 


Declination (J2000) 


8 


-30°39'40.74" 


11°55'55.55" 


16°06'27.40" 


Proper motion in R.A. 


ii a (masyr -1 ) 


_ 


1.34 ±0.01 


-3.27 ±0.35 


Proper motion in Decl. 


H S (masyr -1 ) 
I 




-25.05 ±0.02 


-1.04 ±0.42 


Galactic longitude (J2000) 


245?2 


19?8 


50?0 


Galactic latitude (J2000) 


b 


-4?5 


48?3 


2?1 


Proper motion in Gal. long. 


Hi (masyr- 1 ) 




-21.57 d 


-0.60 


Proper motion in Gal. lat. 


fi b (masyr- 1 ) 




-12.80 d 


-3.38 


Distance 


d (kpc) 


~ 0.6 


1.02 ±0.05 


8.3 ± 1.4 


Characteristic age of the ms pulsar 


t c (Myr) 


210 


250 


110 


Spin-down age of the ms pulsar 


T b (Myr) 


100 


210 


80 


Mass of the ms pulsar 


M A (Mq) 


1.34 


1.33 


1.44 


Mass of the companion 


M b (M q ) 


1.25 


1.35 


1.39 


Current semi-major axis 


A cur (Rq) 


1.26 


3.28 


2.80 


Current orbital eccentricity 


&cur 


0.0878 


0.274 


0.617 



a Burgay ot al. 2003; Lyne et al. 2004. 

^Wolszczan 1991; Stairs ct al. 2002; Konacki, Wolszczan, & Stairs 2003; Arzoumanian, Cordes, & Wasserman 1999. 

c Hulse & Taylor 1975; Taylor ot al. 1976; Taylor, Fowler, & McCulloch 1979; Taylor & Weisberg 1982, 1989; Damour & Taylor 1991; Arzoumanian 
ct al. 1999, Wex et al. 2000. 
d Derived quantities. 



the results presented here to derive tighter constraints. 
We also compare our predicted spin tilts with those de- 
rived based on the geometric model of Jenct & Ransom 
(2004) for pulsar B's eclipses and find excellent agree- 
ment. We note that our predicted spin tilts can also be 
used along with polarization observations to anticipate 
the most likely timing of pulsar A's disappearance due 
to geodetic precession as a function of the radial veloc- 
ity (which in the long term will be the only remaining 
unknown property) of the system. 

In addition to PSR J0737-3039 (Section 2), we analyze 
our current understanding of the evolutionary history of 
the other two close currently known DNS systems in the 
Galaxy. For PSRB1534±12, we use the masses, orbital 
characteristics, and measured proper motion to derive 
progenitor and kick constraints, as well as spin tilt pre- 
dictions (Section 3). For PSRB1913±16, we extend the 
analysis of Wex et al. (2000) to take into account the 
possibility that the helium star progenitor of the pulsar 
companion could have been filling its Roche lobe. We 
also derive the probability distributions of the kick mag- 
nitude imparted to the last-born NS (Section 4) . We end 
with a summary of our main conclusions, a discussion of 
the implications of our analysis for the relation of kicks 
to spin orientations, and a comparison of our results with 
other recent studies. 

2. THE RECENT EVOLUTIONARY HISTORY OF 
PSRJ0737-3039 

In this Section, we use the scintillation velocities mea- 
sured by Ransom et al. (2004) to first determine the spa- 
tial velocity of PSR J0737-3039 with respect to a frame 
of reference locked to the Galaxy and next follow the mo- 
tion of the system in the Galaxy backwards in time until 
the time of its formation. We identify crossings of the 
Galactic plane as possible birth sites and, by subtracting 
the local Galactic rotational velocity from the total sys- 
temic velocity at birth, determine the system's post-SN 
peculiar velocity. This velocity is subsequently used to 
tighten the progenitor constraints derived in Paper I. 



The physical parameters of PSR J0737-3039 relevant 
to this investigation are summarized in Table 1. We par- 
ticularly note that we use the estimate Tf, = 100 Myr of 
the time since pulsar A left the maximum spin-up line as 
an upper limit for the age of the system (for details, see 
Paper I). 

2.1. Present Systemic Velocity 

The interstellar scintillation measurements carried out 
by Ransom et al. (2004) yield two of the three systemic 
velocity components required to describe the space mo- 
tion of PSR J0737-3039. Both components are measured 
in the plane perpendicular to the line-of-sight and with 
respect to the Earth. The first one has a magnitude of 
~ 96.6 ± 3.7 km s -1 and is directed along the line con- 
necting the ascending and descending node of the pul- 
sar's orbit, while the second one has a magnitude of 
~ 103.1 ± 7.7kms~ 1 and is directed perpendicular to 
the nodal line. In what follows, we will refer to these two 
components as V"" and V 1 -, respectively. Since the statis- 
tical errors of 3.7 km s -1 and 7.7 km s -1 do not account 
for unknown systematic errors, we here adopt the values 
— 100 km s -1 and V 1 - = 110km s _1 , and assume an 
uncertainty on both components of ±20kms~ 1 (Ransom 
2004, private communication). 

Before we pass on from the scintillation velocity com- 
ponents to velocity components with respect to a frame 
of reference locked to the Galaxy, it is interesting to ex- 
amine the velocity components V a and V$, and the corre- 
sponding proper motion components ji a and [is, in right 
ascension (a) and declination (S). The velocity compo- 
nents V a and V$ are obtained from VH and V 1 - by means 
of the transformation formulae 

V a = yll cosft- V 1 - shift, 1 
V 5 = shift ±V^- cosft, J 

where ft is the unknown longitude of the ascending node 
measured right-handed around the line-of-sight, from the 
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Fig. 1.— Variation of the velocity components of PSR J0737-3039 
in right ascension a (solid lines) and declination <5 (dashed lines) as 
a function of the unknown longitude of the ascending node f2. The 
dark and light gray bands indicate the errors in V a and Vg corre- 
sponding to errors of 10 and 20kms~ 1 in the scintillation velocity 
components and V^. The right-hand axis shows the magni- 
tude of the proper motion components fi a and fig, for a distance 
of 0.6 kpc. 



East through the North. The proper motion components 
fi a = cos 5 da /dt and [is = dS/dt are then given by the 
relations 



V„ 



kms 

Vs 
kms~ 



= 4.74 



masyr - 
masyr" 



(2) 



d 

,kpc y 
4 ' 74 (kpc, 

where d is the distance of the binary from the Sun. 

The variations of V a (or fi a ) and Vg (or fig) are shown 
in Fig. 1 as a function of the unknown $7. The dark and 
light gray bands indicate the errors in V a (or fi a ) and 
Vg (or [is) corresponding to errors of 10 and 20 kms -1 
in VW and V ± . The magnitude of V a and Vs ranges 
from approximately — 175 kms -1 to 175 kms -1 , while, 
for a distance of 0.6 kpc, the magnitude of fi a and fig 
can be anywhere between —60 masyr" 1 and 60 masyr -1 
(see also Ransom et al. 2004). A 100 mas proper motion 
should then be detectable in less than 2 years time. It 
is furthermore clear that, even with the extended error 
bars of ±20 kms -1 , future measurements of the proper 
motion will considerably constrain the possible values of 
the longitude of the ascending node O. 

Since the scintillation velocity measurements only give 
the magnitude of U" and V ± , the motion of the binary 
projected on the plane perpendicular to the line-of-sight 
should, in principal, be considered in four different di- 
rections corresponding to four different combinations of 
U 11 = ±100 km s -1 and V 1 - = ±110kms -1 . This addi- 
tional degree of freedom may however be incorporated in 
the unknown longitude of the ascending node. To illus- 
trate this, we consider two sets of scintillation velocity 

components {V\ ,V±-) and (uj',^) whose orientations 
in the plane perpendicular to the line-of-sight are deter- 
mined by the angles f2i and ^2- From Eqs. (1), it then 
follows that the two sets of scintillation velocity com- 
ponents give rise to equal velocity components in right 



ascension and declination when fl 2 
AO is determined by 



ill + Af2, where 



cos(AQ) 



uj' - v« V 2 A 



+ 



(U 2 ")2 + (U/)2 



v 1 - v" 

sin(Aft) = 1 2 



vi 



(3) 



Hence, the cases (V^V^) = (-100 kms -1 , 110 kms" 1 ), 
(V\V^) = (100 kms -1 , -110 kms -1 ), and (V^V^) = 
(— 100 kms -1 , — HOkms -1 ), can be obtained from the 
case (V^V^) = (100 kms -1 , 110 kms -1 ) by translat- 
ing f2 over angles of 275.5°, 95.5°, and 180°, respectively. 
For the remainder of the paper, we therefore restrict our- 
selves to presenting the results associated with the veloc- 
ity components U" — 100 kms -1 and V 1 - = 110 kms -1 . 

2.2. Kinematic History and Age 

In order to trace the Galactic motion of PSRJ0737- 
3039 back in time until the time of its formation, we first 
need to determine the three velocity components of the 
current systemic velocity in a reference frame locked to 
the Galaxy. For this purpose, we choose a right-handed 
frame OXYZ with the origin at the Galactic center and 
with the XF-plane coinciding with the Galactic plane. 
We adopt the direction from the Sun to the Galactic 
center as the positive direction of the X-axis and the 
direction of the Sun's motion due to the Galactic rotation 
as the positive direction of the Y-axis. 

The components (Vx,Vy,Vz) of the current systemic 
velocity with respect to OXYZ are obtained by the vec- 
tor transformation 




Vq + J ■ 




(4) 



where Vq is the Galactic velocity of the Sun, J = 
d(X,Y,Z)/d(d,a,5) is the Jacobian of the coordinate 
transformation linking (X,Y,Z) to (d,a,S), and V r is 
the unknown radial velocity. The Galactic velocity of 
the Sun in turn is the vector sum of the Sun's circular 
velocity around the Galactic center and the Sun's pecu- 
liar velocity with respect to the local standard of rest: 



Okms -1 
V Q = ( 221 kms -1 
Okms -1 



lOkms -1 
5 kms -1 
7 kms -1 



(5) 



(Binney & Dehnen 1997, Bienayme 1999). It is evident 
from Eqs. (1) and (4) that the kinematic evolution of 
PSR J0737-3039 depends on two unknown parameters: 
fl and V r . 

The variations of Vx and Vz as functions of f2 are 
displayed in Fig. 2, for V r = -1500, -1250, -1000,..., 
1500 kms -1 . It is interesting to note that for V r < 
—250 kms -1 , Vx is always positive so that the motion 
of the system backwards in time is away from the Galac- 
tic center. In addition, for Q < 25° and Q > 300°, V z is 
positive for all considered values of V r so that the motion 
of the system backwards in time is also away from the 
Galactic plane. As we will see below, this behaviour to- 
gether with the current position of the system below the 
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Fig. 2. — Variation of the velocity components Vx and Vz of 
PSR J0737-3039 as a function of the unknown longitude of the 
ascending node Q, for V r -values ranging from -1500 to 1500 km s — 1 
in increments of 250kms — . 



Galactic plane has important consequences for the past 
motion of PSR J0737-3039 in the Galactic potential and 
the identification of possible birth places. 

Given the present position and velocity components of 
PSR J0737-3039 in the adopted reference frame, we in- 
tegrate the equations governing the system's motion in 
the Galaxy backwards in time and look for possible birth 
sites as a function of SI and V r . Since DNS systems are 
thought to form from primordial binaries with initially 
two massive component stars, which are known to have 
a very small vertical scale height of 50-70 pc, we assume 
the primordial binary to be born in the Galactic disk. 
Wex et al. (2000) furthermore noted that both the or- 
bital velocity before the birth of the first NS and the 
center-of-mass velocity imparted to the system due to 
the first SN explosion are typically of the order of only 
a few lOkrns -1 (see also Brandt & Podsiadlowski 1995; 
Pfahl et al. 2002). However, Pfahl et al. (2002) used a 
population synthesis technique to construct a probability 
distribution of post-SN peculiar velocities of high-mass 
X-ray binaries and found a low-probability tail extending 
all the way up to ~ 200kms _1 . The possibility that the 
first SN explosion in PSR J0737-3039 produced a kick 
of a few 100 km s -1 to the binary's center-of-mass can 
therefore not be entirely ruled out. In order to properly 
treat these cases, the effect of both the first and second 
SN explosion on the space motion of the binary must be 
taken into account, which requires the introduction of 
three additional free parameters (magnitude of the first 
SN kick, direction of the first SN kick relative to that of 
the second kick, and time expired between the two kicks). 
In view of this added complexity and the low probability 
of the high systemic velocities after the first SN explo- 
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Fig. 3. — Kinematic age of PSR J0737-3039 as a function of 
the unknown radial velocity V r and the unknown longitude of the 
ascending node f2, for both case A (top panel) and B (bottom 
panel). 



sion, we here leave aside these more complicated and 
rarer cases and assume the present kinematic properties 
of PSR J0737-3039 to be entirely due to the second SN 
explosion. In particular, we assume that the system was 
still close to the Galactic plane at the time of the second 
SN explosion and that its pre-SN peculiar velocity was 
small in comparison to the Galactic rotational velocity 
at the formation site. The post-SN systemic velocity of 
the binary is then given by the sum of the local Galactic 
rotational velocity and the peculiar velocity due to the 
second SN explosion. 

Following Wex et al. (2000), we use the Galactic po- 
tential of Kuijken & Gilmore (1989) to calculate the 
motion of PSR J0737-3039 in the Galaxy backwards in 
time for all possible il-values between 0° and 360°, and 
K-values between -1500 and 1500 km s -1 . Each time 
the binary crosses the Galactic plane (Z=0), we check 
whether the crossing occurred less than 100 Myr in the 
past and whether the crossing point is within 15kpc from 
the Galactic center. If so, the crossing is considered as a 
possible birth site and the post-SN peculiar velocity V pec 
is determined by subtracting the local Galactic rotational 
velocity from the system's total systemic velocity 2 . 

In cases where the binary crosses the Galactic plane 
within 1 kpc from the Galactic center, we assume that the 
pre-SN binary belonged to the disk- like component of the 
Galactic bulge (see, e.g., Kuijken & Rich 2002). Since the 
time required by a massive binary to evolve into a DNS 
is at most of the order of ~ 20 Myr, much smaller than 
the typical age of 'pure' Galactic bulge stars, we exclude 

2 Note that Table 3 in Wex et al. (2000) contained an error in 
the coefficients j3\ for A = 2, 3; where the correct values are j3\ = 
for both A = 2 and 3. We corrected this error and confirmed that 
it does not quantitatively affect the results presented by Wex at al. 
(2000) for PSRB1913+16 (see also Section 4). In particular, the 
resulting relative errors on the times of the Galactic plane crossings 
and the associated peculiar velocities are less than ~ 6%. 
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Fig. 4. — Cross sections of Fig. 3 showing the variation of the kinematic age Tki n of PSR J0737-3039 as a function of V r , for different 
values of Q. The top and bottom panels correspond to case A and B, respectively. 



birth sites in the bulge at Z 7^ 0. Hence, all disk cross- 
ings within 100 Myr in the past and within 15kpc from 
the Galactic center are considered as possible birth sites 
and are assumed to follow the disk-like rotation curve 
determined by the adopted Galactic potential. 

The time in the past at which the system crosses the 
Galactic disk provides an estimate for the time since 
the formation of the DNS. Since some of the Galactic 
trajectories cross the Galactic plane more than once, 
these kinematic ages are not uniquely determined by ft 
and V r . In particular, for Okms -1 < V r < 550km s -1 
and 60° < ft < 280°, we find that the system may 
have crossed the Galactic plane twice within the last 
100 Myr; while for SOkms" 1 < V r < 400km s" 1 and 
240° < ft < 265° it may have even crossed the Galactic 
plane as much as three times. Since the latter situation 
only occurs for just a few, greatly fine-tuned trajectories 
at ages very close to the adopted upper limit of 100 Myr, 
we leave aside these rather rare cases and focus on the 
first and second Galactic plane crossings. We will refer 
to these crossings as case A and case B, respectively. 

The dependence of the kinematic age Tkin on both ft 
and V r is displayed in Figs. 3 and 4: Fig. 3 shows the 
kinematic ages for all considered ft- and V r -values using 
a linear gray scale (color in the electronic edition), while 
Fig. 4 gives a more detailed view of the dependence of 



Thin on V r for some specific values of ft. 

For case A, the majority of the Galactic disk crossings 
occur less than 20 Myr in the past. These young ages are 
associated with longitudes of the ascending node between 
~ 30° and ~ 260°, and radial velocities larger than ap- 
proximately -1000 kms -1 . Galactic disk crossings occur- 
ring more than 20 Myr in the past are typically associated 
with Si-values smaller than ~ 70° or larger than ~ 230°, 
and with 14-values between -200 kins" 1 and 700 kins" 1 . 
The range of possible V^-values becomes narrower with 
increasing kinematic ages and is most restricted when the 
age of the system is close to the adopted upper limit of 
100 Myr. There are furthermore almost no disk crossings 
associated with negative Upvalues when Tkin ^ 20 Myr. 
The origin of this systematic behavior is related to the 
dependence of the X- and Z-componcnts of the current 
peculiar velocity on V r and ft (see Fig. 2): the Galac- 
tic position of the system is such that negative ^.-values 
correspond to orbits where the first disk crossing occurs 
either too far in the past (in excess of even the character- 
istic age of 210 Myr) or too far away from the Galactic 
center (at radii in excess of 15kpc). The disk crossings 
associated with negative Upvalues and Tkin S 20 Myr oc- 
cur due to the suitable orientation of the orbit in space: 
the ft-values (60° < ft < 230°) are such that the Z- 
component of the current peculiar velocity is negative, 
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Fig. 5.— Peculiar velocity of PSR J0737-3039 as a function of 
the unknown radial velocity V r and the unknown longitude of the 
ascending node Q. The top and bottom panels correspond to case A 
and B, respectively. 



so that motion backwards in time is towards the Galac- 
tic plane. In addition, Lorimer et al. (2004) recently 
estimated the age of the system to be most likely be- 
tween 30Myr and 70Myr. If the age is indeed so con- 
strained, the range of possible O-values would be limited 
to fi < 80° and fi > 240°, and the range of possible 
K-values to -200 kms" 1 <V r < 800 kms" 1 . In view of 
the current uncertainties we will perform our analysis for 
the upper limit of 100 Myr, but the above age constraints 
may be kept in mind. 

Case B disk crossings, on the other hand, only occur 
for Vr-values between Okms -1 and 550km s _1 , and re- 
values between 60° and 280°. The associated kinematic 
ages are always longer than 20 Myr, so that constraining 
the age to 30-70 Myr does not impose severe additional 
limits on the ranges of Q- and ^--values. 

In Fig. 5, we show the peculiar velocities V pec asso- 
ciated with the Galactic plane crossings at ages shown 
in Fig. 3. The peculiar velocities range from 90kms _1 
to 1550 kms -1 and from 120 kms -1 to 800 kms -1 , for 
case A and B, respectively. However, when Tkm ^ 
40 Myr, the post-SN peculiar velocities associated with 
case A are also confined to a smaller interval between 
90kms _1 and 800kms _1 . In addition, for a given value 
of £1, the peculiar velocity Vp OC generally increases with 
increasing absolute values of V r . In Section 2.3, we will 
use the post-SN peculiar velocities to further constraint 
the properties of the pre-SN progenitor of PSRJ0737- 
3039 derived in Paper I. 

We conclude this section by noting that, if PSR J0737- 
3039 is indeed born in the Galactic disk, Figs. 3 and 5 
show that the future determination of may yield im- 
portant constraints on its age and post-SN orbital ve- 
locity, even without knowing the system's present radial 
velocity. 



2.3. Progenitor Constraints 

In this Section, we extend the analysis presented in Pa- 
per I by adding the kinematic history of PSR J0737-3039 
to the constraints resulting from the orbital dynamics of 
asymmetric, instantaneous SN explosions. The two main 
differences with paper I are that (i) the problem is now 
expressed as a function of the two unknown parameters 
fi and V r , and (ii) the spin-down age Th — 100 Myr is 
now used as an upper limit for the age instead of as an 
actual age estimate. The first difference implies that the 
constraints will also become functions of £1 and V r , while 
the second difference affects the amount of orbital evo- 
lution that may have taken place since the birth of the 
system. 

As in Paper I, we assume that the post-SN orbit 
evolved solely under the influence of gravitational radi- 
ation and determine the semi-major axis A and orbital 
eccentricity e right after the second SN explosion by nu- 
merically integrating the orbital evolution equations de- 
rived by Junker & Schafer (1992) backwards in time. 
For each combination of Q and V r , the integration is 
terminated at the kinematic age Tki n corresponding to 
a case A or case B Galactic disk crossing. The result- 
ing post-SN orbital separations and orbital eccentricities 
range from A = 1.26 R Q and e = 0.088 (r kin = OMyr) to 
A = 1.54 R Q and e = 0.12 (r fcm = 100 Myr). 

Under the assumption that the pre-SN binary orbit 
is circular (see Paper I), the pre- and post-SN binary 
parameters are related by the conservation laws of orbital 
energy and orbital angular momentum. The relations 
take the form 

v k +V ( f + 2 V k V Q cosd 
= G(Ma + Mb)(-^-j\, (6) 

Al [V% sin 2 9 cos 2 <p + (V k cos 9 + V Q ) 2 
= G(M A + M B )A(l-e 2 ) (7) 

(e.g., Hills 1983; Kalogera 1996), where G denotes the 
gravitational constant, Aq the pre-SN orbital separation, 

M and V = [G (M A + M ) /A Q } 1/2 the helium star's 
mass and pre-SN orbital velocity, and Vk the magnitude 
of the kick velocity imparted to pulsar B at birth. The 
angles 9 and <fi describe the direction of the kick veloc- 
ity imparted to pulsar B: 9 is the polar angle between 
T4 and Vq, while <f> is the corresponding azimuthal an- 
gle measured in the plane perpendicular to Vq so that 
4> = tt/2 coincides with the direction from pulsar A to 
pulsar B. For a given kick- velocity magnitude and a given 
set of pre- and post-SN binary parameters, the polar an- 
gle 9 £ [0, 7r] is uniquely determined by Eq. (6). The 
azimuthal angle <fr £ [0, 2tt], on the other hand, is only 
determined through the square of its cosine so that, for a 
given kick-velocity magnitude, four solutions to Eqs. (6) 
and (7) exist for each set of pre- and post-SN binary 
parameters. 

A necessary set of conditions for the existence of so- 
lutions of Eqs. (6) and (7) is that < sin 2 9 < 1 and 
< cos 2 <p < 1. Together with the requirements that the 
binary must remain bound after the SN explosion, that 
the post-SN orbit must pass through the position of both 
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stars at the time of the SN explosion, and that the SN 
kick must reproduce the post-SN orbital parameters A 
and e, these conditions impose strong constraints on the 
possible pre-SN values of A and M as well as on the 
magnitude of the kick velocity 14 that have already been 
described in detail in Paper I (see also Fryer & Kalogera 
1997; Kalogera & Lorimer 2000). In addition to these 
constraints, the kick velocity and the pre- and post-SN 
binary parameters must now also be compatible with the 
post-SN peculiar velocity derived by following the motion 
of the system in the Galaxy backwards in time. From 
Eqs. (3) and (34) in Kalogera (1996), it follows that the 
post-SN peculiar velocity is related to the kick velocity 
and to the pre- and post-SN binary parameters as 



V z = 

pec 



M B M 



(M A + M )(M A + M E 



G(M - M B )M A 



M Q A 



+ 



G(M q -Mb)(M -2M b )M a 2 
M Mb Aq \ 



(8) 



We require the peculiar velocity obtained from this equa- 
tion to be within the error bars of the peculiar velocity 
obtained by following the motion of the system back- 
wards in time. In order to estimate the error bars on the 
latter, we neglect the contributions resulting from un- 
certainties in the adopted Galactic potential and assume 
them to be entirely due to the errors in the scintillation 
velocity measurements (±20 kms -1 ). 

In Figs. 6-8, we show the ranges of M , A , and V4 
values satisfying all the imposed constraints as functions 
of the unknown radial velocity V r , for case A and case 
B, and for different specific values of the longitude of the 
ascending node ft. The parameter space for case A disk 
crossings is displayed in the upper 9 panels of the figures, 
whereas the parameter space for case B disk crossings is 
displayed in the lower 6 panels. The gray scale (color 
in the electronic edition) indicates the kinematic ages 
associated with different combinations of ft and V r (cf. 
Fig. 3). 

For each value of ft, there is a finite range of Upvalues 
leading to physically acceptable solutions for the pre-SN 
progenitor of PSR J0737-3039. For case A disk crossings, 
the range is widest for ft ~ 160°, where — 1000 km s -1 < 
V r < 1100 kms -1 . The actual range of V r -vahies leading 
to physically acceptable solutions is therefore somewhat 
smaller than the range considered for the Galactic mo- 
tion calculations presented in Section 2.2. For case B 
disk crossings, the 14-range is widest when ft ~ 80° or 
ft ~ 260°, where Okms -1 < V r < 550 km s -1 . Hence, 
for case B, the physical constraints on the progenitor of 
PSR J0737-3039 do not put any further constraints on 
the radial velocity besides those already found from the 
system's kinematic history (see Fig. 3). 

In Paper I, the mass of pulsar B's helium star progen- 
itor was constrained to be between 2.1 M Q and 4.7 M . 
These limits correspond to the lowest mass for which a 
helium star is expected to form a NS and to the high- 
est mass for which mass transfer from pulsar B's helium- 
star progenitor to pulsar A is expected to be dynamically 
stable. Despite the dependence of the analysis presented 
here on the two unknown parameters ft and V r , the range 
of possible Mo-values is much more tightly constrained 
than in Paper I for many combinations of ft and V r . This 
particularly applies to (ft, V r ) pairs for which the kine- 



matic age is close to the age of 100 Myr adopted there 
(see, e.g., the case A solutions associated with ft = 0° 
and 80 Myr < r kin < 100 Myr in Fig. 6). The range of 
available helium star masses M furthermore tends to be 
most constrained for small values of the radial velocity. 

For a given post-SN orbital separation A and eccentric- 
ity e, the possible values of the pre-SN orbital separation 
Aq range from A(l — e) to A(l + e) (Flannery & van den 
Heuvel 1975). In Paper I, we used the values of A and e 
obtained by integrating the equations governing the or- 
bital evolution due to gravitational radiation backwards 
in time up to an age of 100 Myr to constrain Ao to be 
between 1.36 Rq and 1.72 Rq. These small pre-SN or- 
bital separations imply that the helium star progenitor 
of pulsar B is very likely overflowing its Roche lobe (see 
also Dewi & van den Heuvel 2004). Here we take into 
account that the age of PSR J0737-3039 may be shorter 
than 100 Myr. Consequently, the lower limit on Aq may 
actually be even smaller than 1.36 Rq since less orbital 
evolution may have taken place. An absolute lower limit 
on Aq is given by A cur (l — e cur ) = 1.15 -R©, correspond- 
ing to the lower limit associated with the post-SN orbital 
parameters for the minimum age r^in = Myr. Since the 
increase in A and e due to the orbital evolution back- 
wards in time is largest when T^ n = Tb = 100 Myr, the 
value of 1.72 Rq provides an absolute upper limit on Ao. 
Hence, for the pre-SN helium star mass range quoted 
above, the conclusion reached in Paper I that the imme- 
diate progenitor of pulsar B is most likely overflowing its 
Roche lobe remains valid. The range of possible pre-SN 
orbital separations associated with given values of ft and 
V r is furthermore most restricted for the youngest kine- 
matic ages found: the width of the allowed interval is 
approximately 0.25 Rq when Tki n < 20 Myr and 0.36 Rq 
when Tkin — 100 Myr. 

In Paper I we derived a constrained range for the kick- 
velocity magnitude between 60 kms -1 and 1560 km s . 
The lower limit arises from the requirement that a kick 
counteracts the SN mass loss and leaves the post-SN 
system in a tight orbit with a sufficiently small eccen- 
tricity; the upper limit arises from the requirement that 
the binary remains bound. Here we find that, for a 
given radial velocity, the magnitude of the kick veloc- 
ity is much better constrained: typically to an interval 
that is about ~ 500 kms -1 wide. If the parameter space 
is restricted to kinematic ages close to 100 Myr, as in 
Paper I, the range of allowed kick-velocity magnitudes 
furthermore narrows from 60 kms -1 < V r < 1560 kms -1 
to 60 kms -1 < V r < 710 kms -1 for case A disk crossings 
and to 60 kms -1 < V r < 850 kms -1 for case B disk cross- 
ings. In addition, the kick velocity is smallest when V r 
is small and increases with increasing absolute values of 
V r (see also Wex et al. 2000). Since the post-SN peculiar 
velocity Upe C generally increases with increasing absolute 
values of V r (see Fig. 5), this behavior is consistent with 
Eq. (8). The highest kick velocities furthermore tend to 
be associated with the youngest ages. 

For case A, we find an absolute lower limit on the kick 
velocity imparted to pulsar B of ~ 60 kms -1 and an ab- 
solute upper limit of ~ 1660kms -1 ; while for case B 
the absolute lower and upper limits are ~ 60 kms -1 and 
~ 1390 kms -1 . The lower limits are in agreement with 
the lower limit derived in Paper I and with the lower limit 
derived by Dewi & van den Heuvel (2004). This confor- 
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Fig. 6. — Limits on the mass Mo of PSR J0737-3039B's pre-SN helium star progenitor as a function of V T , for various values of f2 which 
were chosen as to provide the clearest and most informative graphical representation. The upper 9 panels are for case A; the lower 6 panels 
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(cf. Fig. 3). 



PULSAR KICKS AND SPIN TILTS IN CLOSE DOUBLE NEUTRON STARS 



9 



el 



ia 



20 



3D 



4D 



50 



70 



AO 



9Q 



1D0 



1, 3 

■f 1.4 

1.2 
1.S 

£ 1.4 

1.2 
1.8 

1.# 

f 1.4 
1.2 
I.B 

1.4 

1.2 
1,3 

2 

1.4 
1.2 



— 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1— 

" Caw M!- P " 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


— ■ — l — i — i — i — ■ — | — i — i — i — 1—| — i — i — i — i — | — i — i — i — i — | — i— 

- Cqhc a, H - 25" 

■ : 

■ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


— i — | — ■ — i — i — i — | — i—i — i — i — | — i — i — i — i — | — i — ■ — i — i — | — i— 

~ Oaee a, ft - 50" 

■ 

i 1 i i i i 1 i i i i | i i i i | i i i i | i 


- Cose a, ft - 30" " 

1 *1 


- Caw *r - ]W ' 


~ Cqw a, £) — 240" 




i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i 


i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i 


i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 1 i 


" Case a, D - TKf 

\ ; 


" Caw ^ D - 295° 

M 


~ taw a, a - 320° 

1 




Cose a, H - 6C/ 1 " 

1 

i 1 i . i i 1 i . i i 1 i . i i 1 i . i i 1 i 


- Case b„ a - sty ' 

; a 

■ i .... i ■ ... i ■ ... i ■ ... i ■ 


' Casa B, fi - 1 40° 

1 

■ i .... i ■ ... i ■ ... i ■ ... i ■ 


" Coua Bj, Q - 17{) n 

. i .... i .... i .... i .... i . 


' Com B, ft - 235" 

: M 

. i .... i .... i .... i .... i . 


~ Casa B, fl - 250 D 

II 

. i .... i .... i .... i .... i . 



-100D -500 MD 100D -10QD -SOO 5D0 10DD -1OD0 -SOD 5QO 1 □□□ 
V r [km/a} V r fkm/n) V r (km/s) 

Fig. 7. — Limits on the pre-SN orbital separation Aq of PSR J0737-3039 as a function of V r , for various values of Q. The upper 9 panels 
are for case A; the lower 6 panels for case B. The gray scale (color in the electronic edition) indicates the kinematic ages associated with 
different combinations of O and V r (cf. Fig. 3). 
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Fig. 8. — Limits on the magnitude Vfe of the kick velocity imparted to PSR J0737-3039B at its formation as a function of V r , for various 
values of Q. The upper 9 panels are for case A; the lower 6 panels for case B. The gray scale (color in the electronic edition) indicates the 
kinematic ages associated with different combinations of fl and V r (cf. Fig. 3). 
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Fig. 9. — Minimum value of the polar angle 9 between the pre-SN 
orbital velocity and the kick velocity imparted to PSR J0737-3039B 
at its formation as a function of V r (solid lines, bottom axis) and Q 
(dashed lines, top axis). Case A and B are represented by black and 
light gray lines (red and blue in the electronic edition), respectively. 
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Fig. 10. — Correlations between the mass M of PSRJ0737- 
3039B's pre-SN helium star progenitor and the kick velocity V fc 
imparted to its center-of-mass at the time of its formation, for case 
A and Q = 160° . The different gray shades (colors in the electronic 
edition) correspond to different values of Vr- 



mity arises because of the moderate post-SN peculiar ve- 
locities associated with the Galactic disk crossings found 
for \V r \ < 100 km s -1 (see Fig. 5), which impose no addi- 
tional constraints on the system in comparison to those 
already imposed in Paper I. Higher absolute values of the 
radial velocity would impose additional constraints and 
shift the minimum kick velocity to higher values. Ran- 
som et al. (2004) furthermore derived a minimum kick 
velocity of 105 km s -1 by assuming that the scintillation 
velocity component V of the current systemic velocity 
is entirely due to the kick imparted to the second-born 
NS (i.e., that the post-SN peculiar velocity is unaffected 
by the motion in the Galaxy). If we adopt a value of 
V x = 90kms _1 , corresponding to the lowest post-SN 
peculiar velocities found from the Galactic motion calcu- 
lations (see Fig. 5), in conjunction with an upper limit of 
4.7 Mq on the mass of pulsar B's helium star progenitor, 
the procedure adopted by Ransom at al. (2004) yields 
a minimum kick velocity of 75kms _1 . This value is in 
close agreement with the minimum kick velocity derived 
above. 

The upper limits derived for the magnitude of the kick 
velocity are determined by the condition that the binary 
must remain bound after the SN explosion. The condi- 
tion is expressed by the inequality 



V k < 



G{M A 



Mn 



An 



1/2 



1 



M A + M B 
M A + Mo 



1/2' 



(9) 



(e.g. Brandt & Podsiadlowski 1995, Kalogera & Lorimer 
2000). For a given set of NS masses Ma and Ms, the up- 
per limit increases with increasing values of Mq and with 
decreasing values of Aq. Since the absolute lower limit 
on Arj derived here is somewhat lower than in Paper I, 
the absolute upper limit of ~ 1660 km s -1 (case A) on Vk 
is somewhat larger. 

The constraints on the direction of the kick velocity 
imparted to pulsar B at birth are illustrated in Fig. 9, 
where we show the minimum value 9 m in of the angle 9 
between the kick velocity and the helium star's pre-SN 



orbital velocity for which all the imposed constraints may 
be satisfied 3 . The solid lines show the values of 6 m in as 
a function of V r , while the dashed lines show the values 
of Brain as a function of f2. Case A and case B are distin- 
guished by black and light gray lines (red and blue in the 
electronic edition), respectively. The values of 9 min are 
fairly insensitive to the value of O. The dependency on 
V r , on the other hand, is much stronger: 9 m i n reaches a 
minimum value of ~ 115° for |V^-| < 200 km s -1 and then 
increases with increasing absolute values of V r . The ten- 
dency of the kicks to be directed opposite to the orbital 
motion therefore becomes stronger with increasing abso- 
lute values of the radial velocity (and thus with increasing 
values of the kick velocity, see Fig. 8). 

It is also interesting to examine how the constraints 
on the helium-star mass Mq and kick magnitude Vk cor- 
relate. For both cases A and B, all pairs of S7 and V r 
yield results qualitatively similar to the ones illustrated 
in Fig. 10. It is evident that, for a given radial velocity, 
higher masses correlate with lower kicks. This is consis- 
tent with our expectations, since stronger SN mass loss 
requires a smaller kick to achieve the same effect (satisfy- 
ing all of the constraints). Also as expected, as the abso- 
lute value of V r increases, the constrained band shifts to 
higher and higher kick magnitudes. Finally, we note that 
these correlations are derived based on the constraints 
relevant to the formation of PSR J0737-3039 and do not 
necessarily reflect a correlation inherent to the kick mech- 
anism. However, a slope for this unknown correlation 
(because of the physical origin of the kick) similar to the 
one shown in Fig. 10 would require minimal fine-tuning 
of parameters for the formation of PSR J0737-3039, and 
hence would be favored statistically. 

2.4. Isotropic Kick Probability Distribution 

In Paper I, we derived a probability distribution for 
the magnitude of the kick velocity imparted to pul- 

3 The maximum possible value of 8 is usually close to 180° so 
that no additional constraints arise from considering 6 max - 
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Fig. 11. — Probability distribution functions for the magnitude of the kick velocity imparted to PSR J0737-3039B at the time of its 
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sar B, under the assumption that all kick directions are 
equally probable. Here, we adopt the same assumption 
and derive such probability distributions for each com- 
bination of O and V r leading to plausible progenitors of 
PSR J0737-3039. As mentioned before, for either case A 
or B, each pair (fi, V r ) determines the kinematic age and 
hence the post-SN orbital parameters (A, e). Further- 
more, for a given kick- velocity magnitude V/., a pair of 
kick-orientation angles (9, cf>) corresponds to a finite set 
of pre-SN properties (Mq,Aq). Therefore, for each pair 
of (ft, V r ) and for a given value of Vk, we can calculate 
the likelihood of a pair (9, </>) that satisfies the Mq,Aq 
constraints as 



A 



v r ,n( 



)( p\Vk) = A(9,(p\Vk,Vr,n) 
3 sin" 



2vr 2 

= 0, 



if M {6,4>) and A (9,4> 

satisfy the constraints 
if not, 



(10) 



where the polar angle 9 is determined by Eq. (6), and j 
denotes the number of allowed (M ,A ) solutions asso- 
ciated with (9, cf>) for the considered Vk- 

Using Bayes' theorem, the likelihood in Eq. (10) is 
transformed into a probability that, for a given pre-SN 
orbital separation Aq and helium-star mass Mo, the kick- 
velocity imparted to pulsar B had a magnitude Vk- 

Pv rt n{V k \e,<f>) 



_ iV r ,o(Vfc)Ay r ,o(MU fc ) (n) 

(e.g. Bevington & Robinson 2002). Here Pv r ,n(Vk) is 
the probability that, for a given V r and Q, the kick ve- 
locity had a magnitude 14 when no particular set of pre- 
and post-SN binary parameters is imposed. Since, in 
this case, there is no reason to favor any kick-velocity 
magnitude over another, we set Pv r ,n(Vk) — 1. 

The total probability that the kick velocity imparted 
to pulsar B had a magnitude Vk is obtained by deriving 
a distribution function Pv,.,n{Vk\9, <p) for each admissible 
(Aq, Mo) pair and performing the integration 



Pv r ,n( v k\S)= I P^ 
Js 



'v r ,n(V k \9,^)d9d4>. (12) 



Here, S denotes the region of (9, 0)-values corresponding 
to the entire region of admissible Aq- and Mo- values in 
the constrained (^4cb Mo)-parameter space. 

To obtain a set of kick- velocity distributions associated 
with different values of V r , finally, we assume the longi- 
tude of the ascending node to be uniformly distributed 
in the plane perpendicular to the linc-of-sight, and in- 
tegrate the distribution functions Py r , n(Vk\S) over all 
possible values of Q: 

1 f 27r 

Pv r (Vk\S) = — Pv r .n(Vk\S)dn. (13) 
27T Jo 

The probability distribution functions Pv r (Vk\S) are 
displayed in the left-hand panels of Fig. 11 by means 
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of a linear color scale varying from light to dark gray 
with increasing PDF-values (in the electronic edition in- 
creasing PDF values are indicated by yellow, green, or- 
ange, red, and blue). We stress that the figure depicts a 
set of individual one-dimensional kick-velocity distribu- 
tions associated with different Upvalues, and not a single 
two-dimensional distribution function of Vk and V r . The 
right-hand panels of the figure show the probability dis- 
tribution functions (PDFs) associated with some specific 
values of V r . 

For case A disk crossings, the majority of the distribu- 
tions peak at a single value of Vk, which increases with 
increasing absolute values of V r (cf. Fig. 8). This single- 
peaked behaviour is in contrast to the behaviour found 
in Paper I, where the kick-velocity distribution showed 
a primary peak around 150kms _1 , and a smaller sec- 
ondary peak at ~1200-1300kms _1 . The disappearance 
of the double-peaked structure is associated with the ad- 
ditional constraints imposed by the post-SN peculiar ve- 
locities derived from the Galactic motion calculations: 
since kick velocities of the order of 150 km s" 1 give rise 
to very different post-SN peculiar velocities than kick ve- 
locities of the order of 1200-1300 kms -1 , the extra con- 
dition effectively eliminates one of the two peaks in favor 
of the other one. When \V r \ < 400 km s -1 , the peak 
in the distribution functions occurs at kick velocities of 
d00-300kms-\ while for \V r \ > 1000 kms" 1 it occurs 



at kick velocities in excess of 1500 km 



For intermedi- 



ate values of \V r \ (~ 600-800 km s _1 ), a hint of a double- 
peaked structure is still present, although the two peaks 
are much more evenly matched in height and occur much 
closer to each other than those found in paper I. 

For case B disk crossings, the distribution functions al- 
ways show a single-peaked behaviour. The most proba- 
ble kick velocity again increases with increasing absolute 
values of V r and, for V r < 500 km s -1 , ranges from ap- 
proximately 150 km s -1 to 300 km s -1 . For larger abso- 
lute values of the radial velocity, the most probable kick 
velocity is about ~500kms _1 , although only few case B 
disk crossings occur for these high radial velocities. 

Since the probability of a given Vk depends on the con- 
straints derived for Aq and Mq, we tested the robustness 
of the distribution functions to some of the assumptions 
adopted in the previous sections. In particular, we varied 
the constraints on Aq by considering both a lower upper 
limit of 50 Myr and a higher upper limit of 150 Myr on the 
age of the system 4 . We find that neither the lower nor the 
higher age yields significant changes in the distribution 
functions, so that our results are insensitive to reasonable 
changes in the constraints derived for A . In addition, 
the lower limit of 2.1 Mq on the mass Mq of pulsar B's he- 
lium star progenitor was determined by the lowest mass 
for which a non-Roche-lobe-filling helium star is expected 
to form a NS. However, since the pre-SN system is most 
likely undergoing a mass-transfer phase which was initi- 
ated after the cessation of core helium burning, the lower 
limit of 2.1 M should actually be applied to the mass of 
the helium star at the onset of Roche-lobe overflow. The 
mass Mq at the time of the second SN explosion may 
then actually be somewhat lower than 2.1 Mq. We have 

4 Recall that the age affects the range of possible Ao-valucs 
through the constraint A(l — e) < Aq < A(l + e), where A and e 
are the post-SN orbital parameters. 



therefore reconsidered the derivation of the distribution 
functions for a lower limit on Mq of 1.4M Q and did not 
find any significant changes in the PDFs. The reason 
for this is twofold. Firstly, the increase in the admissible 
(A , M )-parameter space resulting from a decrease of 
the lower limit on Mq is small in comparison to the total 
size of the parameter space. A decrease in the minimum 
kick velocity associated with a decrease in the minimum 
helium star mass (cf. Paper I) therefore has a very low 
weight in the construction of the PDFs. Secondly, for 
Mq < 1.8 Mq the minimum kick velocity is no longer 
determined by the lower limit on Mq but by the con- 
straints on the post-SN peculiar velocity Vp 0C . This is 
readily seen from Eq. (8) when one uses the inequalities 
Vpec > 90kms-\ A < 1.72 R Q , and A < 1.54 R Q . Fi- 
nally, the upper limit on Mq was varied by considering 
different critical mass ratios separating dynamically sta- 
ble from dynamically unstable Roche-lobe overflow. The 
behavior of the PDFs again did not show any significant 
changes. 

2.5. Tilt Probability Distributions 

According to our current understanding of the forma- 
tion of DNS binaries, the mass-transfer phase responsible 
for spinning up pulsar A to millisecond periods is also 
expected to align pulsar A's spin-axis with the pre-SN 
orbital angular momentum axis. Depending on its mag- 
nitude and orientation, the SN kick imparted to pulsar B 
at the time of its formation may tilt the post-SN orbit 
with respect to the pre-SN orbital plane, thus causing a 
misalignment between pulsar A's spin axis and the post- 
SN orbital angular momentum axis. The tilt angle A is 
related to the pre- and post-SN binary parameters and 
to the magnitude and direction of the kick velocity as 

1-1/2 



cos A 



A M A + M B 

~A~q M A + Mq 



(1 - e 2 ) 



Vk , 
— cos 
Vq 



1 



(14) 



(Kalogera 2000). The dependence of the tilt angle on the 
direction of the kick velocity is thus entirely through the 
polar angle 9 with no contribution whatsoever from the 
azimuthal angle <f>. 

In order to derive a set of probability distribution func- 
tions Pv r (X\S) for the tilt angle A, we use Eq. (14) to 
eliminate Vk from Eqs. (6) and (7). After the elimina- 
tion, the derivation of the distribution functions is simi- 
lar to the derivation of the the kick-velocity distributions 
P Vr {V k \S) outlined in Section 2.4. The resulting PDFs 
are presented in Fig. 12 in a similar fashion as the kick- 
velocity distributions shown in Fig. 11. 

For case A disk crossings, the tilt-angle distributions 
generally show a single peak which becomes wider and 
moves to larger tilt angles with increasing absolute val- 
ues of V r . For \V r \ < 400 kms -1 , the peak occurs at tilt 
angles A < 30°; while for \V r \ > 1000 km s -1 it occurs 
at A > 135° although it becomes much less pronounced. 
For intermediate Upvalues (~ 600 km s -1 ), a hint of a 
double-peaked structure is present. Tilt angles close to 
A w 90° are furthermore strongly disfavored regardless 
of the value of the radial velocity because they require 
greatly fine-tuned kick magnitudes and directions satis- 
fying Vk cos 9 ~ —Vq (sec Eq. 14). Case B disk cross- 
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Fig. 12. — Probability distribution functions for the misalignment angle A between PSR J0737-3039A's spin axis and the post-SN orbital 
angular momentum axis, for case A and B. The left-hand panels show the entire set of PDFs associated with all admissible V r -values by 
means of a linear color scale which varies from light to dark gray with increasing PDF-values (in the electronic edition increasing PDF 
values are indicated by yellow, green, orange, red, and blue), while the right-hand panels show the PDFs associated with some specific 
W-values. For clarity, the curves in the right-hand panels are offset from each other by an arbitrary amount. 



ings, on the other hand, yield distributions which typ- 
ically peak at tilt angles A < 30°. The peak is most 
pronounced when V r < 400 km s -1 and flattens some- 
what for V r > 400 km s -1 . In the latter case, the dis- 
tributions still favor tilt angles below 70°. As for the 
kick-velocity distributions, these results are very insen- 
sitive to the assumptions adopted in the derivation of 
the A - and Mo-constraints, for both case A and case B. 
The overall behavior of the tilt-probability distributions 
is furthermore in excellent agreement with the analyti- 
cal spin-tilt distributions in coalescing binaries with two 
compact objects derived by Kalogera (2000). 

We conclude this section by noting that the geometri- 
cal model constructed by Jenet & Ransom (2004) to ex- 
plain the flux variations of PSR J0737-3039B predicts tilt 
angles of 16° ± 10° or 164° ± 10°. The PDFs presented in 
Fig. 12 favor the lower tilt angle when \V r \ < 400 kms -1 
and the higher tilt angle when \V r \ > 800 km s -1 . The 
alternative solutions of 82° ± 16° or 98° ± 16° derived 
by Jenet & Ransom (2004) are strongly disfavored by all 
distribution functions. As noted by the authors, these 
solutions are also incompatible with the misalignment 
angle between the pulsar's spin and magnetic dipole axes 
derived by Demorest et al. (2004). 

2.6. Non-isotropic kicks 

The probability distribution functions for the kick- 
velocity magnitude and the spin misalignment angle A 



presented in the previous sections were derived under the 
assumption that all kick directions are equally probable. 
However, in recent years a possible alignment of the SN 
kick with the NS spin axis has been discussed, prompted 
mainly by the alignment of pulsar proper motion vec- 
tors with jet axes in the Crab and the Vela pulsars (Lai, 
ChernofF, & Cordes 2001; Romani 2004). 

Assuming that the NS spin axis is determined by the 
axis of the progenitor helium star and that this axis was 
roughly aligned with the pre-SN orbital angular momen- 
tum axis due to mass transfer, the alignment inferred by 
these observations possibly indicates a preference for kick 
directions perpendicular to the pre-SN orbital plane. In 
order to examine the effect of such preferred kick direc- 
tions, we reconsider the derivation of the kick-velocity 
and tilt-angle distributions under the assumption that 
the kick direction is restricted to be within two oppositely 
directed cones with opening angle £ p and with axes par- 
allel to the pre-SN orbital angular momentum axis (i.e., 
polar kicks). For this purpose, it is convenient to intro- 
duce the angle £ £ [0, ir] between the kick-velocity vector 
and the pre-SN orbital angular momentum vector as a 
measure for the spin-kick alignment. The angle is related 
to the kick-direction angles 9 and <f> by the equation 

cos£ = sin# coscj), (15) 

so that kicks are directed within the above mentioned 
cones with opening angle £ p when 

cos£ p < |sin# cos</>| < 1 (16) 
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Fig. 13. — Limits on the polar angle £ between the pre-SN 
orbital angular momentum axis and the kick velocity imparted to 
PSR J0737-3039B at its formation as a function of V r (solid lines, 
bottom axis) and Q (dashed lines, top axis). Values of £ smaller 
than 90° are lower limits (£ m i n ); values of £ larger than 90° are 
upper limits (^ raoI ). Case A and B are represented by black and 
light gray lines (red and blue in the electronic edition), respectively. 



(Kalogcra 2000). We also note that since, for a 
given kick- velocity magnitude Vk, four solutions (#,</>) 
of Eqs. (6) and (7) exist for each set of admissible pre- 
and post-SN binary parameters, there is also a two-fold 
degeneracy in the solutions for £. 

The range of £ angles available for the formation of 
PSRJ0737-3039 is shown in Fig. 13 as a function of V r 
and f2. Curves below 90° correspond to lower limits, 
£ m i„, and curves above 90° to upper limits, £, m ax- Case A 
and B are indicated by black and light gray lines (red and 
blue in the electronic edition), respectively. It is inter- 
esting to note that it would be impossible to understand 
the formation of PSR J0737-3039, if the kick were very 
closely aligned or anti-aligned to the pre-SN orbital an- 
gular momentum (that could also determine the newborn 
NS spin). Instead a minimum deviation from the pre-SN 
orbital angular momentum axis of 25°-30° appears to be 
necessary to satisfy all the imposed constraints. 

In order to illustrate the effects of polar kicks, we con- 
sidered the case where £ p = 30°. The reduction in the 
available parameter space associated with the restriction 
of the kick direction yields viable solutions for case A 
only when \V r \ < 300 kms -1 , and for case B only when 
Okms -1 < V r < 500 kms -1 . In addition, the mass of 
pulsar B's helium star progenitor becomes constrained 
to Mo ~ 2.1-2.5 Mq, the pre-SN orbital separation to 
Aq « 1.1-1.5 Rq, and the magnitude of the kick velocity 
to V k « 100-600 kins" 1 . 

In Figs. 14 and 15, we show the kick-velocity and tilt- 
angle distribution functions in the case of polar kicks 
with £ p = 30°. Both sets of distributions clearly show a 
much more pronounced single peak than when the kicks 
are assumed to be distributed isotropically (cf. Figs. 11 
and 12). Depending on the radial velocity, the most prob- 
able kick velocity ranges from approximately 200kms _1 
to 550 km s -1 , and the most probable tilt angle from 15° 
to 45°. This range of tilt angles is compatible with the 
tilt angle of 16° ± 10° predicted by Jenet & Ransom 



(2004). 

3. THE RECENT EVOLUTIONARY HISTORY OF 
PSRB1534+12 

Unlike PSR J0737-3039, the relativistic binary ra- 
dio pulsar PSRB1534+12 has an accurately measured 
proper motion with a known direction in right ascen- 
sion and declination (see Table 1), so that its kinematic 
history and progenitor constraints depend only on the 
unknown radial velocity V r . In order to derive these con- 
straints, we adopt the spin-down age t\, = 210 Myr as an 
upper limit to the age of the system. The other physical 
parameters of PSRB1534+12 relevant to the derivation 
are summarised in Table 1. 

Following the same arguments as for PSR J0737-3039 
in Section 2.2, we assume that PSRB1534+12 was close 
to the Galactic plane at the time of the second SN ex- 
plosion and that its pre-SN systemic velocity was almost 
entirely due to Galactic rotation. The possible birth sites 
of the DNS are then obtained by tracing the Galactic mo- 
tion of the system backwards in time as a function of the 
unknown radial velocity V r . We find that, within the im- 
posed age limit of 210 Myr, no Galactic plane crossings 
occur when V r < — 280 km s" 1 , whereas up to four cross- 
ings may occur when V r > — 280 km s -1 . Since four disk 
crossings only occur for relatively few and rather fine- 
tuned Galactic trajectories, we leave these cases aside 
and focus on the possible birth sites associated with the 
first (case A), second (case B), and third (case C) Galac- 
tic plane crossings. 

The kinematic ages and post-SN peculiar velocities as- 
sociated with the Galactic plane crossings are shown in 
Fig. 16 as functions of the radial velocity V r . Case A 
gives rise to a wide range of ages between ~ 1 Myr and 
210Myr, and peculiar velocities between 80kms _1 and 
1500 km s -1 . Cases B and C, on the other hand, yield 
kinematic ages of at least ~ 90 Myr and ~ 130 Myr; and 
post-SN peculiar velocities of 100-220 km s _1 and 80- 
160 km s -1 , respectively. 

The post-SN orbital separation A and orbital eccentric- 
ity e at the times of Galactic disk crossings are obtained 
by numerical integration backwards in time of the equa- 
tions governing the evolution of the orbit under the influ- 
ence of gravitational radiation. The resulting post-SN or- 
bital parameters range from A = 3.28 Rq and e = 0.274 
when T fcin = OMyr to A = 3.36 R Q and e = 0.282 
when Thin — 210 Myr. These post-SN orbital parame- 
ters together with the post-SN peculiar velocities impose 
constraints on the pre-SN progenitor of PSRB1534+12. 
These are derived in a similar way as for the progenitor 
of PSR J0737-3037 (except that the problem depends on 
only one free parameter, V r ). The results of our analysis 
are summarized in Fig. 17. 

Contrary to PSR J0737-3039, which was most likely in 
a state of mass transfer just before the second SN explo- 
sion, PSRB1534+12 may have been detached as well as 
semi-detached at the time the second NS was born. In 
order to illustrate this, the solutions for which no mass 
transfer takes place at the time of the second SN explo- 
sion are indicated by the dark gray (red in the electronic 
edition) regions in Fig. 17, while the light gray (blue in 
the electronic edition) regions indicates the additional 
solutions that become accessible when the possibility of 
mass transfer is taken into account. For the latter solu- 
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Fig. 14. — Probability distribution functions for the magnitude Vk of the kick velocity imparted to PSR J0737-3039B at the time of 
its formation, for polar kicks restricted to two oppositely directed cones with an opening angle of £ p = 30° and with axes parallel to the 
pre-SN orbital angular momentum axis (cf. Fig. 11 for more details). 



tions we adopt the same assumption as before that the 
system may survive the mass transfer phase and form a 
DNS only if the mass ratio Mq/Ma is smaller than 3.5 
(e.g., Ivanova et al. 2003). 

For case A disk crossings, detached pre-SN binary con- 
figurations exist only for V r > 250 kms . For smaller 
radial velocities, the system is always undergoing mass 
transfer from the progenitor of the second-born NS to 
the first-born NS. The pre-SN mass of the helium star 
forming the second NS is constrained within 2.1 Mq - 
8Mq. The lower limit again corresponds to the lowest 
mass for which a helium star is expected to form a NS in- 
stead of a white dwarf, while the upper limit corresponds 
to the highest mass for which a helium star is expected 
to form a NS rather than a black hole (see, e.g., Fig. 1 
in Belczynski, Kalogera, & Bulik 2002; and Table 16.4 
in Tauris & van den Heuvel 2004). The divide between 
the dark and light gray (red and blue in the electronic 
edition) regions at 4.7 Mq corresponds to the adopted 
critical mass ratio Mq /Ma = 3.5 for dynamically stable 
mass transfer. The allowed mass range of pre-SN helium 
star masses is most constrained for \V r \ < 200 kms -1 
when 2.1 Mq < Mo < 4M©. Lower and upper limits on 
the pre-SN orbital separation are given by 2.4 i? Q and 
4.3 Rq. 

The behavior of the kick-velocity magnitude is simi- 
lar as for PSR J0737-3039: the kick velocity increases 
with increasing absolute values of V r and, for a given 
radial velocity, it is generally constrained to an interval 
that is less than ~ 600 kms -1 wide. When all possible 



pre-SN binary configurations are considered, the kick ve- 
locity magnitudes range from 15kms to 1350 kms . 
When configurations for which the helium star is over- 
flowing its Roche lobe are excluded, the range narrows 
to 230 kms -1 < Vjt < 1350 kms" 1 , so that the minimum 
kick velocity is higher than when the possibility of Roche- 
lobe overflow is taken into account. The larger minimum 
kick velocity is required to compensate the larger effect of 
the mass lost from the system during the SN explosion: 
since the minimum Mq for these systems is 4.7 Mq, at 
least 57% of the total pre-SN mass is lost from the sys- 
tem. The minimum and maximum kick velocities are 
furthermore in good agreement with those derived by 
Fryer & Kalogera (1997) and Dewi & Pols (2003). 

The polar angle 6 between the kick velocity and the 
helium star's pre-SN orbital velocity is always larger than 
95°, so that the kicks tend to be directed opposite to 
the orbital motion. When the solutions are restricted to 
those where no Roche-lobe overflow is going on at the 
time of the second SN explosion, the direction of the 
kick velocity is constrained to 9 > 140°. These values 
are again in good agreement with those derived by Dewi 
& Pols (2003). The angle £ between the kick velocity 
and the pre-SN orbital angular momentum axis, on the 
other hand, takes values between 5° and 175°. When the 
parameter space is restricted to the solutions where the 
binary is detached just before the second SN explosion, 
the range of possible £- values narrows to 55° < £ < 125°. 

For case B disk crossings, all allowed pre-SN binary 
configurations imply that the helium star progenitor of 
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Fig. 15. — Probability distribution functions for the misalignment angle A between PSR J0737-3039A's spin axis and the post-SN orbital 
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Fig. 16. — Variations of the kinematic age (left-hand axis) and 
post-SN peculiar velocity (right-hand axis) of PSRB1534+12 as a 
function of the unknown radial velocity V r , for cases A, B, and C. 
The kinematic ages are represented by black solid lines (red in the 
electronic edition) and the peculiar velocities by light gray lines 
(blue in the electronic edition). 



the second-born NS is overflowing its Roche lobe at the 
time of its SN explosion. The mass of the helium star is 
constrained to 2.1 M© < M < 4.1 M©, the pre-SN or- 
bital separation to 2.4 i?© < Ao < 4.3 R@, and the kick- 
velocity mag nitudetoSOkms- 1 <V k < 390 km s" 1 . The 



angle between the kick direction and the helium star's 
pre-SN orbital velocity is restricted 100° < 6 < 180°, 
and the angle between the kick direction and the pre-SN 
orbital angular momentum axis to 10° < £ < 170°. 

Case C disk crossings, finally, also yield no detached 
solutions for the progenitor of PSRB1534+12. The 
constraints for this case are similar to those found for 
case B, except that the mass Mo is always smaller than 
~ 3.1 Mq and the kick velocity T4 is always smaller than 
270kms~ 1 . 

From these constraints, it follows that, if the case for 
the alignment of kicks with the spin of the NS's progeni- 
tor gets more support in the future, the direct progenitor 
of PSRB1534+12's companion must be overflowing its 
Roche lobe, regardless of which disk crossing corresponds 
to the birth place. In addition, as for PSR J0737-3039, 
the formation of PSRB1534+12 cannot be explained if 
the kick imparted to PSRB1534+12's companion was 
perfectly aligned or anti-aligned with the spin of its pro- 
genitor. 

The constraints derived above may again be used to 
derive probability distribution functions for the magni- 
tude T4 of the kick velocity imparted to the second-born 
NS and for the misalignment angle A between the spin- 
axis of the first-born NS and the post-SN orbital angular 
momentum axis (cf. Sections 2.4 and 2.5). The result- 
ing PDFs are presented in Figs. 18 and 19 as those for 
PSRJ0737-3039 in Figs. 11 and 12. For radial veloci- 
ties of only a few 100 kms -1 , the kick-velocity distribu- 
tions show two closely spaced and fairly evenly matched 



18 



Willems, Kalogera, & Henninger 



+.0 

? 3.5 
I 

2.5 
1KB 



E 

■■ J 




ieo 



130 







Hi 


a 


D 


ii 








\jjk 


□ 


1 j 


am 



ieo 



120 



&o 



-250 D 25D 5DD "TM 1000 -200 -100 D 1DD 



-1DD -3D D 5D 1DD IM 



Fig. 17. — Limits on the pre-SN progenitor of PSRB1534+12 and on the kick velocity imparted to the last-born NS. The left-hand 
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Fig. 19. — Probability distribution functions for the misalignment angle A between PSR B1534+12's spin axis and the post-SN orbital 
angular momentum axis, for cases A, B, and C (cf. Fig. 12 for more details). 



peaks between Vk — 50kms _1 and Vk — 250 km s -1 . 
For higher radial velocities, relevant only to case A, the 
peak(s) shift to larger kick velocities up to a maximum of 
~ 1350 km s -1 . The tilt-angle distributions, on the other 
hand, favor misalignment angles below 30° when \V r \ < 
200 km s" 1 and above 100° when \V r \ > 600 km s" 1 . For 
case A disk crossings, tilt angles close to A « 90° are fur- 
thermore strongly disfavored regardless of the value of 
the radial velocity. For case B disk crossings, tilt angles 
with non- vanishing probabilities are always smaller than 
40°-50°. 

4. THE RECENT EVOLUTIONARY HISTORY OF 
PSRB1913+16 

Besides the knowledge of the proper motion direction, 
PSRB1913+16 has the advantage that the misalignment 
angle between the pulsar's spin axis and the pre-SN or- 
bital angular momentum axis has been determined to be 



around ~ 20°, corresponding to prograde rotation, or 
~ 160°, corresponding to retrograde rotation (Kramer 
1998; Weisberg & Taylor 2002). Wex et al. (2000) used 
this information to derive constraints on the mass of the 
second-born NS's direct progenitor, on the pre-SN orbital 
separation, and on the magnitude and direction of the 
kick velocity imparted to the second-born NS at birth. In 
agreement with the then available numerical simulations 
of rapidly accreting neutron stars, the authors assumed 
that mass transfer from a helium star companion would 
cause the NS to collapse into a black hole and therefore 
excluded Roche-lobe overflowing helium stars as viable 
progenitors of the second-born NS. However, more re- 
cent calculations by Dewi et al. (2002), Ivanova et al. 
(2003), and Dewi & Pols (2003) show that NS binaries 
may survive a helium-star mass-transfer phase, provided 
that the ratio of the helium star's mass to the neutron 
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star's mass is not too extreme (Mq/Ma < 3.5). In this 
section, we therefore revise the constraints derived by 
Wex et al. (2000) in the light of this new information. 
The adopted physical parameters are listed in Table 1. 
The constraints on the tilt angle, for which we consider 
the values A = 18° ± 6° and A = 162° ± 6°, are imposed 
by means of Eq. (14). 

We look for possible birth sites of PSRB1913+16 by 
tracing the system's motion in the Galaxy backwards in 
time up to a maximum age of 80Myr 5 , as a function of 
the unknown radial velocity V r . In agreement with Wex 
et al. (2000), we find that the system may have crossed 
the Galactic plane up to two times. The first Galactic 
plane crossing (case A) occurs at very young kinematic 
ages of ~ 2-4 Myr and gives rise to peculiar velocities 
in excess of ~ 300 km s -1 . The corresponding post-SN 
orbital parameters are A w 2.8 Rq and e « 0.618. The 
second Galactic plane crossing (case B), on the other 
hand, takes place at least ~ 55 Myr in the past and yields 
peculiar velocities of ~ 230-440 km s _1 . The associated 
post-SN orbital separations and eccentricities range from 
A = 3.1 Rq and e = 0.646 to A = 3.2 Rq and e = 0.658. 

The constraints on the pre-SN parameter space acces- 
sible to the progenitor of PSRB1913+16 are shown in 
Fig. 20 as functions of the unknown radial velocity V r . 
As before, the dark gray (red in the electronic edition) 
regions correspond to the solutions for which no mass 
transfer takes place at the time of the second SN explo- 
sion, while the light gray (blue in the electronic edition) 
regions indicate the additional solutions found when the 
possibility of mass transfer is taken into account. The 
constraints for detached pre-SN binary configurations are 
in good agreement with the constraints derived by Wex 
et al. (2000). It is clear however, that when the possibil- 
ity of mass transfer is taken into account, the available 
parameter space becomes much less constrained. 

For case A disk crossings, the range of radial veloci- 
ties for which physically acceptable solutions exist is re- 
stricted to \V r \ < 500-600 km s" 1 when A = 18°, and 
to lOOkms- 1 < \V r \ < ISOOkms- 1 when A = 162°. 
The mass of the second-born NS's direct progenitor is 
constrained to the interval between 2.1 M© and 8M Q , 
and the pre-SN orbital separation to the interval be- 
tween 1.1 Rq and 5.3 Rq for both considered values of 
the tilt angle A. The magnitude of the kick velocity varies 
from ~ 190 km s -1 to ~ 600kms _1 when A = 18° and 
from ~ 580 km s" 1 to ~ 2000 km s" 1 when A = 162°. 
When the solutions are restricted to detached pre-SN 
binary configurations the limits become 300 km s -1 < 
U fc < 600 km s" 1 and 700 km s' 1 < V k < 1250 km s" 1 
for A = 18° and A = 162°, respectively. The direction of 
the kick velocity is most constrained when A = 162°: for 
A = 18° we have 95° < 9 < 175° and 5° < |90°-£| < 85°, 
while for A = 162° we find 160° < 9 < 175° and 
5° < 90° - £| < 15°. When the parameter space is 
restricted to detached pre-SN binary configurations, the 
limits become 145° < 9 < 175° and 5° < |90° - £| < 15° 
for A = 18°, and 165° < 9 < 175° and |90° - £| « 5° for 

5 As for PSRJ0737-3039 and PSRB1534+12, we use the spin- 
down age r;, = 80 Myr instead of the characteristic age r c = 
110 Myr as an upper limit for the age of the system. The max- 
imum amount of orbital evolution that may have taken place since 
the formation of the DNS is therefore somewhat smaller than in 
Wex et al. (2000). 



A = 162°. 

For case B disk crossings, the mass Mq of the second- 
born NS's direct progenitor and the pre-SN orbital sep- 
aration A are constrained to the ranges of values given 
by 2.1 Mq < M < 8M Q and 1.1 R Q < A < 5.3 i?© 
when A = 18°, and to 2.1 M < M < 5M Q and 
2.9 Rq < A < 5.3 Rq when A = 162°. The magni- 
tude of the kick velocity varies between ~ 190 km s _1 and 
~ 530 kms -1 when A = 18°, and between ~ 550 kms -1 
and ~ 850 km s -1 when A = 162°. When the solutions 
are restricted to those where no Roche-lobe overflow oc- 
curs at the time of the helium star's SN explosion, the 
minimum kick velocity associated with A = 18° increases 
slightly to approximately 280 km s -1 , while the range of 
admissible kick velocities associated with A = 162° be- 
comes very tightly constrained to Vk ~ 640-680 km s _1 . 
The direction of the kick is again most tightly constrained 
when A = 162°: for A = 18° we have 80° < 6 < 175° and 
|90°-£| > 5°, while for A = 162° we find 160° < 9 < 175° 
and 5° < |90° - £| < 10°. When the parameter space 
is restricted to detached pre-SN binary configurations, 
9 and £ are constrained to 140° < 9 < 175° and 
5° < 90° - £| < 15° for A = 18°, and to 6 w 170°- 
175° and 90° - £| w 5° for A = 162°. 

Hence, for both case A and B and for both possible 
tilt angles A, the formation of PSRB1913+16 can only 
be explained if the kicks imparted to the second-born NS 
are "elevated" out of the orbital plane by at least 5° . In 
addition, the constraints on 9 and £ imply that a large 
tilt angle of 162° cannot be obtained if the kicks are per- 
pendicular to the pre-SN orbital plane (i.e. aligned or 
anti-aligned with the pre-SN orbital angular momentum 
vector). Instead, for an unambiguous determination of 
a large tilt angle A, our results would be more compat- 
ible with the findings of Deshpande, Ramachandran, & 
Radhakrishnan (1999) and Birkel & Toldra (1997) who 
argue against a spin-kick correlation. 

The inclusion of pre-SN binary configurations for 
which stable mass transfer occurs from a helium star 
to the first-born NS furthermore considerably increases 
the available parameter space for the formation of 
PSRB1913+16. This is especially true for the prograde 
case, A = 18°, for which the kick direction is no longer 
required to be close to the pre-SN orbital plane as was 
found by Wex et al. (2000). The conclusion by Wex et 
al. (2000) that the retrograde case, A = 162°, requires 
more fine-tuning in the kick direction than the prograde 
case remains valid even when the possibility of pre-SN 
Roche- lobe overflow is taken into account. 

The probability distribution functions for the magni- 
tude of the kick velocity imparted to the second-born NS 
in PSRB1913+16 are shown in Fig 21 for both A = 18° 
and A = 162°. The distributions generally show a sin- 
gle peak at kick velocities which increase with increasing 
absolute values of V r . For A = 18°, the most probable 
kick velocity ranges from ~ 300 km s -1 to ~ 600 km s -1 , 
while for A = 162° it ranges from ~ 600 km s" 1 to almost 
~ 2000 km s" 1 . 

As a test, we also constructed tilt-probability distribu- 
tions for PSRB1913+16 by omitting the tilt-angle con- 
straints from the derivation of the admissible parameter 
space. For case A, the distribution functions peak at 
A < 30° when \V r \ < 400 kms" 1 and at A > 150° when 
\V r \ > 900 kms -1 . Intermediate absolute values of the 
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Fig. 20. — Limits on the pre-SN progenitor of PSRB1913+16 and on the kick velocity imparted to the last-born NS. Dark gray (red 
in the electronic edition) regions correspond to solutions for which the pre-SN binary is detached, while light gray (blue in the electronic 
edition) regions indicate the additional solutions associated with mass-transferring systems. 



radial velocity yield predictions for the tilt angle which 
arc incompatible with the observationally derived values. 
For case B, the distributions peak at tilt angles A < 30° 
and become vanishingly small for A > 30° over the en- 
tire range of possible ^.-values (see Fig. 20). A large 
tilt angle of A = 162° can therefore only be explained 



if the system has crossed the Galactic plane only once, 
which would imply that its age is of the order of only a 
few Myr. It is also interesting to note that a tilt angle 
of 18° would be most favored for radial velocities around 
-300 km s- 1 and +400 km s"\ and a tilt angle of 162° 
for any radial velocity smaller than -900kms _1 or larger 
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Fig. 21. — Probability distribution functions of the magnitude of the kick velocity imparted to PSR B1913+16's companion at the time 
of its formation, for both A = 18° and A = 162° (cf. Fig. 11 for more details). In the right-hand panels, the solid lines correspond to 
A = 18° and the dashed lines to A = 162°. 



than +900 km s _1 . 

5. CONCLUSIONS AND DISCUSSION 

In this paper, we have extended the investigation by 
Willems & Kalogera (2004) and have derived a more elab- 
orate and tighter set of constraints for the progenitor of 
PSR J0737-3039 just before the SN explosion that forms 
the second NS as well as for the magnitude and the di- 
rection of the kick imparted to this NS at birth. The 
additional constraints stem from the kinematic history 
of the system which is derived by tracing the system's 
motion in the Galaxy backwards in time and by assum- 
ing that its peculiar velocity results entirely from the 
second SN explosion. In order to trace the system's mo- 
tion backwards in time, we use the scintillation veloc- 
ity measurements carried out by Ransom et al. (2004), 
which yield two velocity components in the plane per- 
pendicular to the line-of-sight. Since the orientation of 
the velocity components is unknown, their orientation 
O in the plane perpendicular to the line-of-sight intro- 
duces the first free parameter in the problem. In order to 
fully describe the system's space motion, the scintillation 
velocity components must furthermore be supplemented 
with the system's unknown radial velocity V r , the second 
free parameter in the problem. 

We find that, depending on the values of and V r , the 
system may have crossed the Galactic plane up to two 
times. We identify these crossings as possible birth sites 
of PSR J0737-3039 and use the times of the crossings as 
estimates for the system's age. It follows that there is 



a wide range of fi- and ^-values for which PSRJ0737- 
3039 may be remarkably young (< 20 Myr). If the system 
has crossed the Galactic plane twice, it is at least 20 Myr 
old regardless of the values of f2 and V r . Our kinematical 
analysis furthermore shows that the post-SN peculiar ve- 
locity must have been at least on the order of ~ 90 kms -1 
and may even have been as high as ~ 1200 kms -1 (see 
Fig. 5). This is in contrast to the work done by Piran 
& Shaviv (2004) who used the present proximity of the 
system to the Galactic plane to derive a statistical 95% 
confidence upper limit of ~ 150 kms -1 on the post-SN 
peculiar velocity. 

Next, we used the age estimates and peculiar veloci- 
ties associated with the Galactic plane crossings to im- 
pose additional constraints on the progenitor and forma- 
tion of PSR J0737-3039 besides those already discussed in 
Willems & Kalogera (2004). We find that despite the two 
degrees of freedom, the pre-SN and kick parameters as 
well as the age become much more constrained for many 
combinations of fl and V r (see Figs. 6-8). An overview of 
the constraints on the mass Mq of pulsar B's helium star 
progenitor, the pre-SN orbital separation Aq, and the 
magnitude Vk of the kick velocity imparted to pulsar B 
at birth as a function of V r for all possible values of f2 is 
given in Fig. 22. The absolute limits covering the whole 
range of admissible Q- and U r -values are summarized in 
Table 2. 

We also examined constraints on the kick direction in 
terms of the angle 9 between the kick direction and the 
pre-SN orbital velocity of pulsar B's helium star pro- 
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TABLE 2 

Summary of the constraints derived for PSR J0737-3039, PSRB1534+12, and PSRB1913+16. For PSRB1534+12, and 
PSRB1913+16 numbers between parentheses correspond to the constraints associated with detached pre-SN progenitors. 



Parameter 


PSRJ0737-3039 


PSRB1534+12 




PSRB19134 


-16 


A f = 18° ± 6° 






A f = 162° ± 6° 


M (M Q r 


2.1-4.7 


2.1-8.0 (4.7-8.0) 


2.1-8.0 (4.5-8.0) 






2.1-8.0 (4.5-8.0) 


A (R Q ) h 


1.2-1.7 


2.4-4.3 (2.5-4.3) 


1.1-5.3 (2.5-5.3) 






1.1-5.3 (2.5-5.3) 


V k (kms- 1 ) c 


60-1660 


15-1350 (230-1350) 


190-600 (300-600) 






580-2000 (700-1250) 


8 (degrees) d 


115-180 


95-180 (140-180) 


95-175 (145-175) 






160-175 (165-175) 


£ (degrees) 


25-155 


5-175 (55-125) 


5-85 or 95-175 (75-85 or 


95-105) 


75- 


-85 or 95-105 («85 or «95) 



Prc-SN mass of the last-born NS's helium star progenitor. 
^Pre-SN orbital separation. 
c Kick- velocity magnitude, 
d 



Angle between the kick direction and the pre-SN orbital velocity of the last-born NS's helium star progenitor. 
Angle between the kick direction and the pre-SN orbital angular momentum axis. 

Misalignment angle between PSRB1913+16's spin axis and the post-SN orbital angular momentum axis. 



j^H ^Hf ^| ' 



500 - 




-1000 -500 500 

V r (km/s) 
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V r (km/s) 
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Fig. 22. — Overview of the constraints on the mass Mo of PSR J0737-3039B's helium star progenitor, the pre-SN orbital separation Aq, 
and the magnitude Vk of the kick velocity imparted to pulsar B at birth, as a function of V r for all possible values of fl. 



genitor, as well as the angle £ between the kick direction 
and the pre-SN orbital angular momentum vector (which 
could be associated with pulsar B's spin axis). We find 
that 9 is always larger than ~ 115°, while £ ranges from 
~ 25° to ~ 155°. The kick imparted to pulsar B at birth 
must therefore have been directed opposite to the orbital 
motion and cannot have been too closely aligned with the 
pre-SN orbital angular momentum axis. 

For each combination of Q and V r yielding viable pro- 
genitors for PSR J0737-3039, we derived probability dis- 
tributions for the magnitude of the kick velocity under 
the assumption that all kick directions are equally proba- 
ble (Fig. 11). The distribution functions exhibit a single 
peak at kick velocities which increase from ~ 100 km s -1 
to ~ 1600-1700 kms -1 with increasing absolute values 
of the radial velocity. This is in contrast to Willems 
& Kalogera (2004) who found a single most probable 
kick velocity of ~ 150 kms . The dependence of the 
position of the peak on the radial velocity found here 
stems from the relation between the kick magnitude and 
the post-SN peculiar velocity which tends to select those 
kick velocities which are able to explain the system's mo- 
tion in the Galaxy (high space velocities cannot be ex- 
plained by low kick velocities and vice versa). However, 
for a given radial velocity, the range of kick velocities 
with non-vanishingly small probabilities found here is al- 



ways much more constrained than in Willems & Kalogera 
(2004). 

Similarly, we derived distribution functions for the mis- 
alignment angle A between pulsar A's spin axis and the 
post-SN orbital angular momentum axis (Fig. 12). Tilt 
angles lower than 30°-50° are clearly favored in the cases 
where (i) the system has crossed the Galactic plane twice 
in the past and (ii) the system crossed the Galactic 
plane once and has a current radial velocity of less than 
~ 500 kms -1 in absolute value. Higher tilt angles must 
be associated with one disk crossing and much higher 
present-day radial velocities of ~ 500-1200 km s -1 in ab- 
solute value. Tilt angles close to 90° are furthermore 
strongly disfavored for any radial velocity V r . Our prob- 
ability distributions are thus compatible with the tilt- 
angle predictions of 16° ± 10° and 164° ± 10° made by 
Jenet & Ransom (2004) and with the spin- magnetic 
dipole misalignment angle derived by Demorest et al. 
(2004) which rules out Jenet & Ransom's alternative so- 
lutions of 82° ± 16° or 98° ± 16°. 

In view of recent claims of alignment of NS kicks with 
NS spin axes (e.g. Romani 2004), we also considered 
kick-velocity and tilt-angle distributions for polar kicks 
where the kick direction is confined within two oppo- 
sitely directed cones with an opening angle of 30° and 
with axes parallel to the pre-SN orbital angular momen- 
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turn axis (i.e. £ < 30°). The confinement of the kick 
directions strongly reduces the range of radial velocities 
for which viable progenitors for PSR J0737-3039 may be 
found. Correspondingly, the range of most probable kick 
velocities shrinks to ~ 200-550 km s -1 , and the range of 
most probable tilt angles to 15°-45° (see Figs. 14 and 15). 

Progenitor constraints and isotropic kick-velocity 
and spin-tilt distributions were also derived for 
PSRB1534+12 and PSRB1913+16. Both systems have 
a measured proper motion with a known direction in the 
plane of the sky, so that the derived constraints and prob- 
ability distributions only depend on the unknown radial 
velocity V r . For PSRB1913+16, the knowledge of the tilt 
angle A (Kramer 1998; Weisberg & Taylor 2002) further- 
more puts additional constraints on the possible progen- 
itor systems, which are not yet available for PSR J0737- 
3039 and PSRB1534+12. 

The absolute limits on the pre-SN binary and kick 
parameters of the two systems are summarized in Ta- 
ble 2 (the dependency of the constraints on V r is shown 
Figs. 17 and 20). Contrary to PSR J0737-3039, the pro- 
genitors of these two systems were not necessarily un- 
dergoing a mass-transfer phase at the time the second 
NS was born. The helium star masses Mq, pre-SN or- 
bital separations A , and kick velocities T4 associated 
with detached progenitors are, however, always much 
more constrained than those associated with the full set 
of detached and semi-detached progenitors. In particu- 
lar, the possibility of pre-SN Roche-lobe overflow con- 
siderably relaxes the constraints on the progenitor of 
PSRB1913+16 derived by Wex et al. (2000). 

The constraints on the direction of the kick im- 
parted to the second-born NS in PSRB1534+12 and 
PSR B 19 13+ 16 yield very similar conclusions as for 
PSR J0737-3039: the kicks are generally directed op- 
posite to the orbital motion and cannot have been too 
closely aligned with the pre-SN orbital angular momen- 
tum axis. 

The kick-velocity and spin-tilt distributions for 
PSRB1534+12 and PSRB1913+16 (Figs. 18, 19 and 21) 
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also show qualitatively the same behavior as for 
PSRJ0737-3039. The most probable kick velocity de- 
pends on the unknown radial velocity and, depending 
on the value of V r , can be anywhere between the de- 
rived lower and upper limits for the kick imparted to the 
second-born NS. Tilt angles smaller than 30°-50° are 
again favored when (i) the system crossed the Galac- 
tic plane more than once and (ii) the system crossed 
the Galactic plane once and has a current radial ve- 
locity smaller than ~ 250 km s -1 in absolute value for 
PSRB1534+12 and smaller than ~ 500 km s" 1 in abso- 
lute value for PSRB1913+16. 

For conclusion, we point out that the set of observa- 
tional measurements and constraints available for dou- 
ble neutron stars systems are sufficient to greatly limit 
the properties of the progenitors. As a follow-up inves- 
tigation it will be therefore interesting to compare these 
strongly constrained progenitor properties to synthesis 
studies of binary populations and double neutron star 
formation. 
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